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PRODUCTION  RATES  FOR  ELECTRON  BEAMS  AND  SWARMS  IN  NITROGEN 


1.  Introduction 

The  deposition  of  electron  energy  in  nitrogen  is  of  interest  in  many 
different  fields,  e.g.,  electron  beam  propagation  in  the  atmosphere, 
stopping  of  auroral  electrons,  electron-beam-pumped  lasers.  To  solve  this 
problem  the  cross  sections  for  all  relevant  electron-atom  and  electron- 
molecule  collision  processes  must  be  known.  These  cross  sections  are  then 
incorporated  into  the  Boltzmann  equation,  from  which  the  electron 
distribution  function  is  obtained.  Once  the  distribution  function  is 
found,  production  rates  for  various  processes  are  calculated.  This  report 
provides  additional  details  which  were  not  included  in  a  previous 
publication*^  (denoted  I)  on  the  application  of  the  deposition  model'^"  to 
nitrogen. 

The  form  of  the  Boltzmann  equation  needed  to  solve  the  energy 

deposition  problem  and  the  nitrogen  cross  sections  are  summarized  in  Secs. 

2  and  3.  Section  4  contains  theoretical  results  for  electron  swarms  in 

nitrogen.  For  applied  fields  ranging  from  1  -  300  Td  (1  Townsend  > 

-17  2 

10  ev-cm  ),  the  drift  velocity,  characteristic  energy,  and  production 
rates  for  the  various  processes  are  given.  The  fraction  of  energy  going 
into  different  channels  is  shown.  Results  are  compared  with  those 
obtained  by  assuming  a  Maxwellian  electron  distribution  function.  Section 
5  gives  beam  deposition  results  where  the  beam  electric  field  is  ignored 
and  the  beam  energy  is  varied.  The  energy  expended  per  electron-ion  pair, 
V,  is  presented,  along  with  production  rates.  Finally,  results  are 
compared  with  those  which  assume  the  beam  electrons  are  completely 
stopped. 


Manufcripi  approved  January  4,  1990. 
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2. 


Boltzmann  Equation 

1-5 

The  secondary  electron  distribution  function  is  given  by 


|f(T..)  -  N  |e  [«jK)  v(t*Ej)  £(T,Ej.t) 


2T+I, 


-  (Tj(T)  v(T)  f(T,t) 


-  E 

i 


J  de  a'[e,  e-Ij-lj  v(e)  f(e,t) 


T+I. 


r  ® 

+  J  de  cJ(C,T)  v(e)  f(c,t)  -  a^(T)  v(T)  f(T,t) 

2T+Ij 

+  S(T,t)  -  ||(T,t)  -  ||(T,t)  ,  (1) 

where  f(T,t)  is  the  secondary  electron  density  per  unit  energy 
-3  -1 

(cm  eV  )  for  electrons  with  kinetic  energy  T  and  speed  v(T),  the 
maximum  secondary  electron  energy,  and  N  the  density  of  the  molecules 
which  is  assumed  constant  for  the  times  of  interest.  Populations  of  the 
excited  states  are  negligible.  Equation  (1)  includes  energy  loss  by 
electrons  to  all  excitations  with  a  cross  section  of  and  an 

excitation  threshold  of  E^,  as  well  as  to  ionization,  where  aj(T)  is  the 
total  ionization  cross  section  of  N2  resulting  in  the  i-th  ionization 
continuum  and  (r^(e,T)  is  the  differential  ionization  cross  section. 

In  Eq.  (1),  S(T,t)  is  the  source  term  for  the  generation  of 
electrons  with  energy  T  by  the  beam  electrons.  It  is  assumed  that  the 
high  energy  primaries  leave  the  volume  of  interest  after  at  most  one 
collision  and 

S(T,t)  =  N  Njj(t)  v(T^)  p  »i(Tb»T)  ,  (2) 


where  N|^(t)  is  the  electron  density  of  a  beam  whose  energy  is  and  the 
summation  is  over  all  ionization  channels. 

The  term  D(T,t)  in  Eq.  (1)  represents  the  flux  of  the  secondary 
electrons  driven  by  the  beam-induced  electric  field  and  is  given  by 


D(T,t) 


2Ne^(EQ/N)^T 


2mc' 


3m(v  /N) 

ID 


3T 


T°*^ll  + 


f(T.t) 


(3) 
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where  v  is  the  electron  collision  frequency  for  momentum  transfer  and 
the  electric  field  strength.  The  elastic  collision  term,  H(T,t),  is 


2mNv 


H(T,t)  = 


m 


M 


f(T,t) 


.5T  -  T 

e 


_  T  tMILlH 
3T 


(4) 


where  M  is  the  mass  of  the  target  and  T  is  the  target  temperature  (in 
eV).  For  the  problems  considered  in  this  paper,  3H/9T  is  small  compared 
to  the  other  terms  in  Eq.  (1). 


3.  Cross  Sections 

The  model  consists  of  one  effective  rotational  cross  section,  eight 
vibrational  excitation  cross  sections,  six  triplet  excitations,  four 
singlet  excitations,  eleven  channels  leading  to  pure  dissociation,  four 
dissociative  ionization  channels,  three  pure  ionization  channels,  inner- 
shell  ionization  and  excitation  cross  sections.  The  momentum  transfer 
cross  section  is  also  Included.  A  discussion  of  these  cross  sections  with 
references  to  their  measurements  is  contained  in  I.  Figure  1  shows  the 
total  cross  sections  for  these  processes.  From  these,  the  loss  function 
is  easily  obtained  and  given  in  Fig.  2. 

4.  Electron  Swarm  Results 

To  simulate  an  electron  swarm,  the  beam  source  in  Eq.  (1)  was  set  to 
zero,  some  seed  electrons  were  assumed,  and  a  constant  electric  field  was 
applied  until  an  equilibrium  was  reached.  If  the  field  is  high  enough, 
equilibrium  is  one  in  which  the  distribution  function  grows  at  a  steady 
exponential  rate  in  time  because  of  ionization.  The  range  of  fields 
tested  was  from  1  -  300  Td.  For  lower  fields  the  results  are  sensitive  to 
the  details  of  the  rotational  cross  sections  and  for  the  higher  fields  the 
two-term  approximation  of  the  Boltzmann  equation  becomes  less  accurate. 

For  the  swarm  results,  the  background  gas  temperature  is  T  -  0.025  eV. 

-3  ® 

Figure  3  gives  the  normalized  (Xf(e)de«  1cm  )  distribution  function 
for  several  values  of  the  reduced  electric  field  strength,  E/N  (where  N  is 
the  density  of  the  neutrals).  Figure  4  shows  the  calculated  drift 
velocity  v^  as  a  function  of  E/N.  Figure  5  gives  the  characteristic 
energy  and  average  energy  of  the  electron  distribution  function. 

To  compare  the  results  with  those  obtained  by  assuming  a  Maxwellian 
distribution,  the  average  electron  energy  is  used  to  provide  a 
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temperature,  T  «  2/3e  ,  for  the  Haxvellian.  Figure  6  plots  the 

ionization  and  pure  dissociation  coefficients.  These  results  compare  veil 
with  experimental  measurements  as  shown  in  I. 

Figure  7  shows  the  average  collision  frequency  for  momentum 
transfer.  This  collision  frequency  was  obtained  in  two  different  ways. 
First,  ■  eE/mVj,  where  is  the  drift  velocity.  This  is  the  proper 
average  to  use  in  the  calculation  of  the  conductivity,  a  «  e^n^/mv^. 
Second,  ■  Jf (c)<Tj|j(e)vdc.  Also  shown  on  the  graph  is  the  average 

momentum  transfer  collision  frequency,  ..  ,  calculated  with  f  given  by  a 
Maxwellian  with  the  same  total  energy. 

Figure  8  gives  the  fraction  of  the  energy  deposited  by  the  field  in 
the  various  processes.  Vibrational  excitation  accounts  for  the  majority 
of  absorption  until  about  100  Td  where  excitation  of  the  triplet  states 
begins  to  dominate. 

3 

Figure  9  gives  the  production  rates  (#/cm  -sec)  for  many  of  the 
species  in  the  model. 

Finally,  Fig.  10  compares  the  distribution  function  and  several  of 
the  production  rates  with  those  assuming  a  Maxwellian  of  the  same  energy. 
The  assumption  of  a  Maxwellian  distribution  tends  to  underestimate  triplet 
production,  but  overestimate  ionization.  Newman  and  DeTemple^  have  shown 
similar  results. 

5.  Electron  Beam  Deposition  Results 

For  these  studies  the  electric  field  is  turned  off  and  the  beam 

source  term  S(T,t),  given  in  Eq.  (2),  is  applied.  For  energies  above  the 

1-3 

lowest  threshold,  the  distribution  function  comes  to  an  equilibrium. 

Figure  11  shows  V,  the  amount  of  energy  expended  by  the  beam  in  order  to 

produce  an  electron-ion  pair,  as  the  beam  energy  is  varied.  V  is  nearly 

independent  of  beam  energy  for  energies  above  a  few  kilovolts.  Also  shown 

in  Fig.  11  are  the  values  of  V  obtained  by  assuming  the  source  term  is  a 

delta  function  at  a  fixed  energy;  this  corresponds  to  a  completely-stopped 

1-4  7 

beam  electron.  This  assumption  is  more  accurate  at  lover  energies.  ’ 

Figure  12  gives  the  normalized  equilibrium  distribution  function  for 
a  1  MeV  beam  and  Fig.  13  shows  the  energy  deposited  per  electron-ion  pair 
in  the  various  channels.  Finally,  Fig.  14  shows  the  production  rates,  in 
terms  of  the  number  of  events  per  electron-ion  pair,  for  many  of  the 
processes  in  the  model. 
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6 .  Summary 

Electron  interaction  with  molecular  nitrogen  has  been  studied  for 
1)  electrons  which  are  accelerated  by  a  constant,  uniform  electric  field 
and  2)  high-energy  beam  electrons.  Details  of  the  electron  energy 
deposition  are  given. 
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Attn:  Dr.  Edward  P.  Lee 

Dr.  Thomas  Fessenden 


Lawrence  Livermore  National  Laboratory 
University  of  California 
Livermore,  California  94550 
Attn:  Dr.  Simon  S.  Yu 

Dr.  Frank  Chambers 

Dr.  James  W.-K.  Hark,  L-477 

Dr.  Villiam  Fawley 

Dr.  Villiam  Barletta 

Dr.  Villiam  Sharp 

Dr.  Daniel  S.  Prono 

Dr.  John  K.  Boyd 

Dr.  John  Clark 

Dr.  George  J.  Caporaso 

Dr.  Donald  Prosnitz 

Dr.  John  Stewart 

Dr.  Y.  P.  Chong 

Major  Kenneth  Dreyer 

Dr.  Hans  Kruger 

Dr.  Thaddeus  J.  Orzechowski 

Dr.  Michael  R.  Teague 

Mr.  John  T.  Veir 

Dr.  James  E.  Leiss 
13013  Chestnut  Oak  Drive 
Gaithersburg,  MO  20878 

Lockheed  Missiles  and  Space  Co. 

3251  Hanover  St. 

Bldg.  205,  Dept  92-20 
Palo  Alto,  CA  94304 
Attn:  Dr.  John  Siambis  ^ 

Los  Alamos  National  Laboratory 
P.O.  Box  1663 
Los  Alamos,  NM  87545 
Attn:  Dr.  L.  Thode 

Dr.  H.  Dogliani,  MS-5000 
Mr.  R.  Carlson,  MS-P940 
Dr.  Carl  Ekdahl,  HS-0410 
Dr.  Joseph  Hack 
Dr.  Melvin  I.  Buchwald 
Dr.  David  C.  Hoir 


Maxwell  Laboratories  Inc. 

8888  Balboa  Avenue 
San  Diego,  CA  92123 
Attn:  Dr.  Ken  Vhitham 

McDonnell  Douglas  Research  Laboratories 
Dept.  223,  Bldg.  33,  Level  45 
Box  516 

St.  Louis,  HO  63166 
■Attn:  Dr.  Carl  Leader 

Dr.  Frank  Bieniosek 
Dr.  John  Honig 
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Mission  Research  Corporation 
1720  Randolph  Road,  5.E. 
Albuquerque,  NH  87106 
Attn:  Dr.  Thomas  Hughes 

Dr.  Lawrence  Wright 
Dr.  Kenneth  Struve 
Dr.  Michael  Hostrom 
Dr.  Dale  Welch 

Mission  Research  Corporation 
P.  0.  Drawer  719 
Santa  Barbara,  California  93102 
Attn:  Dr.  C.  Longaire 
Dr.  N.  Carron 

National  Inst,  of  Standards  &  Tech. 
Gaithersburg,  Maryland  20760 
Attn:  Dr.  Mark  Wilson 

Naval  Postgraduate  School 
Physics  Department  (Code  61) 
Monterey,  CA  93940 
Attn:  Prof.  John  R.  Neighbours 
Prof.  Fred  Buskirk 
Prof.  Kai  Woehler 
Prof.  Xavier  Maruyama 

Naval  Surface  Warfare  Center 
White  Oak  Laboratory 
Code  R-41 

Silver  Spring,  Maryland  20903-5000 
Attn:  Mr.  W.  M.  Hinckley 
Dr.  M.  H.  Cha 
Dr.  H.  S.  Uhm 
Dr.  R.  Fiorlto 
Dr.  K.  T.  Nguyen 
Dr.  R.  Stark 
Dr.  H.  C.  Chen 
Dr.  0.  Rule 
Dr.  Matt  Brown 
Mrs.  Carolyn  Fisher  (G42) 

Dr.  Eugene  E.  Nolting  (H23) 

Office  of  Naval  Research 
800  North  Quincy  Street 
Arlington,  VA  22217 
Attn:  Dr.  C.  W.  Roberson 
Dr.  F.  Saalfeld 

Office  of  Naval  Research  (2  copies) 
Department  of  the  Navy 
Code  01231C 
Arlington,  VA  22217 


Office  of  Under  Secretary  of  Defense 

Research  and  Engineering 

Room  3E1034 

The  Pentagon 

Washington,  DC  20301 

Attn:  Dr.  John  MacCallum 

Physics  International,  Inc. 

2700  Merced  Street 
San  Leandro,  CA.  94577 
Attn:  Dr.  E.  Goldman 

Dr.  James  Benford 
Dr.  George  B.  Frazier 
Mr.  Ralph  Genuario 

Princeton  University 
Plasma  Physics  Laboratory 
Princeton,  NJ  08540 
Attn:  Dr.  Francis  Perkins,  Jr. 

Pulse  Sciences,  Inc. 

600  McCormack  Street 
San  Leandro,  CA  94577 
Attn:  Dr.  Sidney  Putnam 

Pulse  Sciences,  Inc. 

2001  Wilshire  Boulevard 
Suite  600 

Santa  Monica,  CA  90403 
Attn:  Dr.  John  R.  Bayless 
Dr.  R.  Adler 

The  Rand  Corporation 
2100  M  Street,  NW 
Washington,  DC  20037 
Attn:  Dr.  Nikita  Wells 
Nr.  Simon  Kassel 

Sandia  National  Laboratory 
Albuquerque,  NM  87115 
Attn:  Dr.  David  Hasti/1272 
Dr.  Collins  Clark 
Dr.  John  Freeman/1241 
Dr.  Charles  Frost 
Dr.  George  Kamin/1274 
Dr.  Gordon  T.  Leifeste 
Dr.  Gerald  N.  Hays 
Dr.  Michael  G.  Nazarakis/1272 
Dr.  John  Vagner/1241 
Dr.  Ron  Lipinski/1274 
Dr.  James  Poukey 
Dr.  Hilton  J.  Clauser/1261 
Or.  Kenneth  R.  Prestwich/1240 
Or.  Kevin  O'Brien 
Or.  Isaac  R.  Shokair 
Or.  J.  Pace  VanOevender/1200 
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Science  Applications  Inti.  Corp. 

5150  El  Camino  Road 
Los  Altos,  CA  94022 
Attn:  Dr.  R.  R.  Johnston 
Or.  Leon  Feinstein 
Dr.  Douglas  Keeley 
Or.  E.  Roland  Parkinson 

Science  Applications  Inti.  Corp. 

1710  Goodridge  Drive 
HcLean,  VA  22102 
Attn:  Mr.  V.  Chadsey 
Dr.  A  Drobot 
Or.  K.  Papadopoulos 
Or.  William  V.  Rienstra 
Dr.  Alan  J.  Toepfer 
Dr.  Alfred  Mondelli 
Dr.  D.  Chernin 
Dr.  R.  Tsang 

Science  Research  Laboratory,  Inc. 

1600  Wilson  Boulevard 
Suite  1200 
Arlington,  VA  22209 
Attn:  Dr.  Joseph  Mangano 
Or.  Daniel  Birx 

Commander 

Space  &  Naval  Warfare  Systems  Command 
PMW-145 

Washington,  DC  20363-5100 
Attn:  CAPT  J.  D.  Fontana 
LT  Fritchie 

SRI  International 
PSO-15 

Molecular  Physics  Laboratory 
333  Ravensvood  Avenue 
Menlo  Park,  CA  94025 
Attn:  Dr.  Donald  Eckstrom 

Dr.  Kenneth  R.  Stalder 

Strategic  Defense  Initiative  Org. 

SDIO/T/DEO 

The  Pentagon 

Washington,  DC  20009-7100 
Attn*.  Lt  Col  R.  L.  Gullickson 
Dr.  D.  Duston 

Titan/Spectron,  Inc. 

P.  0.  Box  4399 
Albuquerque,  NM  87196 
Attn:  Dr.  R.  Bruce  Miller 
Dr.  John  Smith 


Titan  Systems,  Inc. 

2685  Marine  Way 
Suite  1408 

Mountain  View,  CA  94043 
Attn:  Dr.  Kenneth  W.  Billman 

Titan  Systems,  Inc. 

9191  Towne  Centre  Dr. -Suite  500 
San  Diego,  CA  92122 
Attn:  Dr.  R.  M.  Dove 

University  of  California 
Physics  Department 
Irvine,  CA  92664 
Attn:  Dr.  Gregory  Benford 
Dr.  Norman  Rostoker 

University  of  California 

San  Diego,  CA  92110 

Attn:  Dr.  Marshall  N.  Rosenbluth 

University  of  Maryland 
Physics  Department 
College  Park,  MD  20742 
Attn:  Dr.  Y.  C.  Lee 
Dr.  C.  Grebogi 
Dr.  W.  Destler 
Dr.  C.  Striffler 

University  of  Michigan 
Dept,  of  Nuclear  Engineering 
Ann  Arbor,  MI  48109 
Attn:  Prof.  Terry  Kammash  - 
Prof.  R.  Gilgenbach 

Director  of  Research 

U.S.  Naval  Academy 

Annapolis,  MD  21402  (2  copies) 
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Naval  Research  Laboratory 
Washington.  DC  20375-5000 
Code  1220 

Naval  Research  Laboratory 
Washington.  DC  20375-5000 
Code  2630 
Timothy  Calderuood 
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